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Abstract 

Gradual numbers have been introduced recently as a means of extending stan- 
dard interval computation methods to fuzzy intervals. The literature treats mono- 
tonic functions of fuzzy intervals. In this paper, we combine the concepts of grad- 
ual numbers and optimization, which allows for the evaluation of any differentiable 
function on fuzzy intervals, with no monotonicity requirement. 



1 Introduction 

The evaluation of a function on a fuzzy interval has long been accomplished by evalu- 
ating the function on each a-cut of the the fuzzy interval. Fortin, Dubois, and Fargier 
recently introduced the concept of gradual numbers GJQ]]. A fuzzy interval can be 
represented as an interval of gradual numbers, which allows the extension of standard 
interval computation methods to fuzzy intervals. 

Previously published results have extended the practice of analyzing monotonic 
functions on intervals by evaluating the interval extreme points [2|. In the present pa- 
per, we employ optimization techniques on gradual numbers to evaluate differentiable 
functions on fuzzy intervals with no monotonicity requirements. 



2 Classical Interval Computation 

Classical interval analysis can be viewed as seeking upper and lower bounds for the 
value of a operator on interval-valued arguments. For example, when the operator is a 
real-valued function, /: R" h- ► R and x\ £ \x[ ,x^\,...,x n € [x~,x+], the goal is tofindz~ 
andz + such that z=f{x\ ,---,x„) lies in the interval for every x G X = Xi[xf ,Xj ]. 

Warmus (9J, Sunaga O, and Moore |6| first proposed what Lodwick O calls an 
axiomatic approach to interval arithmetic, which defines the following rules for basic 
interval arithmetic: 
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Forx £ [x ,x + ] andy £ [y ,y + ]: 

1. x + y = [x~ + y~,x++y+], 

2. x-y = [x~ -y + ,x + -y~], 

3. ix)> = [min{x~ x y~ ,x~ x y + ,x + x y~,x + x y+},max{x~ x y~,x~ x y + ,x+ x 
y-,x+xy+}], 

4. x-^y = [min{x~ -^-y~,x~ -H-y+,x+ + y~,x + -H-y+},max{x~ -^y~,x~ -H-y+,x+ 

y-,* + -^y + }].0^[y-,y + ]. 

Axiomatic interval arithmetic assumes that all variables are independent, including 
duplicate copies of the same variable. This can lead to overestimations of the bounds 
of functions. Consider, for example, 

f(x) = x — x, x£ [0,1]. 

Using the axiomatic approach yields the result 

[1,0] -[1,0] = [0-1,1-0] = [-1,1]. 

If, however, we simplify the function before applying axiomatic interval analy- 
sis, we get f(x) = x — x = 0, which indicates that [z~,z + ] = [—1,1] is a gross over- 
estimation. 

Simplifying functions to reduce multiple copies of a variable can eliminate the 
overestimation, but is not always possible. For example, there is no way to express 
f(x) = xe* with only one copy of the variable x. 

When the function in question is monotonic, we can appeal to a well-known propo- 
sition to find exact bounds for [z~,z + ] without overestimation. Before stating this re- 
sult, we provide some preliminary definitions. A configuration, co, is an n-tuple which 
lies in X = Xj[xf ,x^]. The set of extreme configurations, 9i , are n-tuples in which ev- 
ery element takes on the value of an interval endpoint. Of course, there are 2" elements 
in H . Now consider the following proposition, which indicates what is known as the 
vertex method of interval analysis in the fuzzy set community. 

Proposition 1: Letx = (x\, ...,x n ) be a tuple of n real interval-valued variables such 
that Xj £ [xj ,x~l\, and let z = f(x\,...,x n ) = [z~,z + ]. If / is continuous and locally 
monotonic with respect to each argument, then 

z~ = min/(<») 

mew 

and 

z + = max/(co). 

This follows from the fact that X = Xi[xf ,x^] is a closed subset of R n , and from 
the Weierstrass theorem that a continuous function attains its maximum and minimum 
on a compact set. The vertex method, while effective, is limited in scope. We are still 
left with over-estimations or underestimations on non-monotonic functions. Take, for 
example, 
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f(x)=x(l-x) for xE [0,1]. (1) 

If we consider the two instantiations of x to be independent, !tf contains the four ele- 
ments {(0,0), (0,1), (1,0), (1,1)}, and the result is 

z- = /(to)- = min /(co) = /(l, 1) = /(0, 0) = /(0, 1 ) = 0, 

mew 

and 

z + =/(co)+ = max/(co) =/(l,0) = 1. 

CfleW 

Inspection reveals that there is no single value of x in the interval [0, 1] for which 
f(x) = x(x — 1) has a value of 1, so [z~,z + ] = [0, 1] is an unsatisfactory solution. We 
might, on the other hand, consider the two instantiations of x to be the same variable. 
Then 9i contains only the two elements {(0,0), (1, 1)}, and the result is 



and 



z- =/(©)- = min /(co) =/(l,l) =/(0,0) =0, 
o>eW 



z + = /(co)+ = max/(co) =/(!,!)= /(0,0) = 0. 



But [z~,z + ] = [0,0] is also wrong because zero is not the only possible value of /. 
Consider, for example, f{\)~j^- We require a broader tool than the vertex method for 
evaluation of non-monotonic functions of interval-valued variables with dependencies. 



3 Optimization and Interval Analysis 

Observe that the values we seek in interval analysis, z~ = mm x£X f(x) and z + — 
max xGX f(x) are the answers to the optimization problems: 

Minimize f(x) (2) 

subject to xi > xj , for i = l,...,n (3) 

— Xi > — x\ , for i = 1, ...,n, 

and 

Maximize f(x) (4) 

subject to Xi > xj , for i = 1 , ...,« (5) 

—Xi > — Xj~, for i = 1, ...,n. 

Lodwick proposed this idea in the nineties with his constrained interval analysis 
(15)), and more recently Kreinovich (|3|) has applied it to technical problems with 
partial information. 

Viewing the interval analysis problem as a mathematical programming problem 
puts at our disposal the rich optimization theory of the last century, which leads to the 
following proposition. 



3 



Proposition 2: Let x — (x\,...,x n ) be a tuple of n variables such that x; € fa ,xf \ 
and let z — f{x\ , ...,x„) such that / is continuously differentiable with respect to each 
X(, Now consider the set of KKT configurations, %7 and %f , to be configurations 
which satisfy the Karush-Kuhn-Tucker conditions for ©-(O and Q respectively. 
Then 

z~ = min /(«) 

<aex~ 

and 

z + = max /(co). 

coesc+ 

Proof: Finding the minimum function value for/(xi , ...,x„) where x 61 = x,[x ( ^,x^] 
can be stated as a mathematical programming problem: 

Minimize z '■ f{x) (6) 
subject to c(i) : x ( — > 0, forj'=l,...,n (7) 
c(n + /) : xt — x,- > 0, for i = l,...,n, (8) 

where x is a 1 x n vector. 

Since constraints c(z') and c(n + /) are lower and upper bounds forx,-, at most one at a 
time of c(z') and c(n + i) will be active. We observe that our mathematical programming 
problem satisfies the linear independence constraint qualifications (i.e. the gradients of 
the active constraints are linearly independent), which means that the KKT conditions 
are necessary conditions for optimality. The details of Karush-Kuhn-Tucker conditions 
are found in any standard optimization textbook, for example Q. 

Now any optimal point will satisfy the following: 



\J f(x)-Xi{xj- Xi ) - ... 

X 




(9) 


X n (x n -x~) + X n+ i(xi -x\) + ... 




(10) 


— + X2n{xb-X+) = 




(11) 


Xi>0, i=l, 


. . . , 2n 


(12) 


Xi-x^>0, z*=l, 


...,« 


(13) 


-x,+x+>0, i=l, 


...,« 


(14) 


Xi{xi-xj)=0, z'=l, 


...,n 


(15) 


X„ + i{—xi+xf) =0, ;' = 1 , 


...,n. 


(16) 



The set of points satisfying (fm >-(fT~5T> are exactly the points in . Note that X = 
x,-[x^,x^] is a closed subset of R", so is a compact set. The Weierstrass theorem 
guarantees that a continuous function attains a maximum and minimum on a compact 
set, so we know /(x) has a minimum on X, which has to be min ;ce „-/(x). 

To find the maximum, we replace f(x) in (O with — /(jc), and find the KKT points, 
which are exactly %f . 
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3.1 Examples of Interval Analysis via Optimization 

For the sake of illustration, recall example (Q]): 

f(x)=x(l-x),xe [0,1]. (17) 

The Karush-Kuhn-Tucker conditions for (JTTJ are 

-2x+l -Xi + X 2 = (18) 

h > 0,7 = 1,2 (19) 

x > 0, (20) 

1 — jc > 0, (21) 

Xix = (22) 

X 2 (l-x) = 0. (23) 

%7 consists of the three points which satisfy (fT8l-([T6]i: 

• x = 0; X = (1,0); f(x)=0 

• x = l; X=(0,l); f(x)=0 

• x= i; X = (0,0); /(*) = £. 

Now zT = min vG5C - f(x) = 0. In this example, it happens that {x £ %^ } = {x G 

SO ^ + = maX xG5C+/( X ) = V 

Optimization yields the desired result, [z~,z + ] = [0, j\ without overestimation or 
underestimation. It should be noted that this accuracy is not without its cost. The vertex 
method tells us that, for monotonic functions / : R" i— » R, function extrema will always 
occur at interval endpoints, so \?(\ = 2 n . In contrast, the number of KKT configura- 
tions of an interval constrained mathematical programming problems will sometimes 
be fewer than the number of extreme configurations, and will occasionally be many, 
many more. In the worst case, the number of points satisfying the KKT conditions will 
be infinite. 



4 Gradual Numbers and Fuzzy Interval Analysis 

Up to this point, we have considered intervals of real numbers, which provide a Boolean 
model for uncertainty. Numbers inside the interval are entirely possible, while numbers 
outside the interval are impossible. It is often the case, however, that the sets we wish 
to model are not so crisply defined. The need frequently arises to move gradually from 
impossible to possible, with values on the outskirts of the interval being somewhat 
possible. Fuzzy set theory accomplishes this by making the bounds of the interval 
softer, so that the transition from an impossible value to an entirely possible value is 
gradual. 

A fuzzy interval, W can be defined by a membership function, /j^ : R i— > (0,1], 
where fj$, (x) indicates to what extent x is a member of W. A standard fuzzy interval is 
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the L — R fuzzy interval, which is defined by a membership function /j^, and reference 
functions Land/?. L: [0,inf) — > [0,1] and/?: [0,inf) — > [0,1] are upper semi-continuous 
and strictly decreasing in the range (0,1], and fjy, is defined as follows: 

{1 forx 6 [w - , w + ], 

M 1 ^) forxov- (24) 

Fortin, Dubois, and Fargier (2) recently proposed that a fuzzy interval can be rep- 
resented as a set of gradual numbers that lie between two gradual number endpoints 
in the same way that a real interval is a set of real numbers that lie between two real 
endpoints. 

Definition 1 : A gradual real number (or gradual number), r, is defined by an as- 
signment function Af : (0, 1] R. It can be understood as a real value parametrized by 
a € (0, 1]. Then for each a, a real value r a is given by A?(a). Let G be the set of all 
gradual real numbers. 

We can describe a fuzzy interval W by an ordered pair of two gradual numbers, 
[vv~,m> + ]. Since a fuzzy interval is a fuzzy set with an upper-semicontinuous mem- 
bership function, we require that A lP - be increasing, that A$+ be decreasing, and that 
5= Ai>+(1)- F° r our results, we will additionally assume that A tI - and A lT ,+ are 
continuous and that A$- (0) and A lP + (0) are defined. 

Before considering arithmetic operations on intervals of gradual numbers, we need 
a definition of arithmetic operations on gradual numbers themselves. If r and s are two 
gradual numbers, then f + s is defined by its membership function ftf+s = A? + 
Subtraction, multiplication, division, minimum and maximum are similarly defined. In 
addition, r is said to be less than s if Sl?((X) < J4s(a) for all a € [0, 1]. Greater than and 
equality are similarly defined. We would like to note that the set of gradual numbers is 
not a fully ordered set. 

Beyond arithmetic operations on gradual numbers, we will find it useful to extend a 
real-valued function f[x\ , ...,x„) : R" i— > R to a gradual function, f(x\,...,x n ) :G" >—> G. 
We define the assignment function of the value of / as follows: 

*-f[ h = /(ri(a)-,r„(a)), (25) 

Va e [0, 1] (26) 

A fuzzy interval W = [w~ ,w + ] is the set of all gradual numbers, w, for which 
w~ < w < w + . If W is an L-R fuzzy interval, then w~ and w + are defined by the 
following assignment functions: 

-%-(«) = w~-L~ 1 (a)a w (27) 
^+(a) = w + +R- 1 (a)$ Wl (28) 

where (respectively R~ l ) is the inverse of L (respectively R) in the part of its 
domain where it is positive. 
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4.1 Earlier Results on Applying Interval Analysis to Fuzzy Inter- 
vals 

The representation of a fuzzy interval as a crisp interval of gradual numbers allows 
us to use tools previously developed for real interval analysis in the analysis of fuzzy 
intervals. 

Fortin, Dubois and Fargier extend the vertex method of interval analysis to intervals 
of gradual numbers Q. They define a fuzzy configuration, £2, which is a tuple of 
gradual numbers, (n,...,r„) as follows: 

Let fuzzy intervals x\ . . .x„ be described by their gradual endpoints, [x J" , x ^ } . . . [x~ , x „ ] . 
A configuration, £2, is an n-tuple of gradual numbers which lies in X = Xi[x[,Xf~]. The 
set of extreme fuzzy configurations, H, are n-tuples in which every element takes on 
the value of an interval endpoint. A fuzzy configuration or extreme configuration eval- 
uated at a is defined as the tuple of each of its gradual components evaluated at a. 
That is D. (a) = {f\ (a) , f„ (a)) . We now turn to the vertex method for fuzzy interval 
analysis. 

Proposition 3: Let I = (x\,...,x n ) be a tuple of n fuzzy interval-valued variables 
such that xt e [x ; ~ , jcj 1 "] . Let / : G" i— » G be the gradual extension of / : R" R as defined 
in (l25l >. and define the interval to be the smallest interval of gradual numbers 

that contains £ = {/(Jci, ...,x„)} for every x G X = Xi[x[ ,x^~]. If / is continuous and 
locally monotonic with respect to each argument, then 

and 

Z+ =max ne ^/(£2). 
The reader is referred to [2| for the proof. 

An example of the vertex method for intervals of gradual numbers will be instruc- 
tive. Let x = [x~ ,x + ] be a fuzzy interval defined by the following gradual endpoints: 

JT(a) = ia, ae[0,l] (29) 

x+{a) = l--a, ae [0,1], (30) 

(31) 

as in Figure Q] 

Let / = x 2 for x&X. To apply the vertex method, we determine that there are two 
elements in H: i~(a) = ioc, a £ [0, 1] andl + (a) = 1 - ia, a G [0, 1]. Now 

/(JT(a)) = V,ae[0,l] (32) 
./(x+(a)) = (1-^oc) 2 (33) 
= 1-a + ^a 2 , ae [0,1]. (34) 
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It is clear that f(x~(a)) < f(x + (a)) for all a e [0, 1], so z~ = f(x") and z + = 
f{x + ). The fuzzy interval z — [z~,z + ] is shown in Figure[2] 

The vertex method, while effective, is limited in scope. We are still left with over- 
estimations, underestimations, or meaningless results on non-monotonic functions. For 
example, letX be defined as above, and let /be the gradual extension of f(x) = x(l —x) 
for x £ X. To apply the vertex method, we determine that there are two elements in H: 
x~(a) = jU, a £ [0, 1] andi + (a) = 1 - \a, a £ [0,1]. Now 

f{x-{a)) = I<x(l-i<x), ae[0,l] (35) 

= I a -ia 2 , ae[0,l] (36) 
= f(x + (a)). (37) 

Then i~=z + , so the fuzzy interval is not an interval at all, but a single gradual number 
defined by the assignment function -%(oc) = jOC— ^(ot) 2 , as in Figure[3] 



5 Main Results 

Fortin, Dubois, and Fargier have used the concept of gradual numbers to extend interval 
techniques to the evaluation of monotonic functions of fuzzy intervals. In this section 
we expand the use of gradual numbers. We employ optimization theory in conjunction 
with gradual numbers to evaluate continuously differentiable non-monotonic functions 
of fuzzy-interval valued variables. 
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5.1 Optimization on Fuzzy Intervals 

Let / : G" t— > G be the gradual extension of a continuously differentiable function / : 
R" ^Ras defined in ((25), and define the interval [z~ ,z + ] to be the smallest interval of 
gradual number that contains z = {f(xi, ...,x n )} for every x G X — Xi[x[ ,Xi~]. 

Observe that the values we seek, zT 
answers to the optimization problems 



and 



Minimize 
subject to 



Maximize 
subject to 



m 

Xi 
Xi 

m 

Xi 
Xi 



> x { , for i = 1 , . 
for i — 1 , . 



<xj, 



<5c1 



for i = 1 , . 
for / = 1 , . 



(38) 
(39) 



(40) 
(41) 



We now define the set of KKT configurations, KlJ and %t, to be configurations 
which satisfy the Karush-Kuhn-Tucker conditions for (|38l-([39l and d40b-(T4Tb respec- 
tively. 

Explicitly, %f is the set of every configurations G X which satisfies the following: 



g(x)-Xi(xi-x[)-... 

■Xn(x„-X~)+K+l(Xi -xf) + ... 

— +hn(Xi-X%) = 

Xi>0, i=l, 

Xi — x^ > 0, i= 1, 

-I,+i+>0, z=l, 

A,;(i,-xr) = 0, z = 1 , 

X„ + /(-I i +x+) = 0, i=l, 



2n 

n 

n 

n 

n, 



(42) 
(43) 
(44) 
(45) 
(46) 
(47) 
(48) 
(49) 



where ,g is the gradual extension of \7/(xi , ...,x„). We define SC + (oc) = {i(a)|lG SC + . 
% ~ (a) is similarly defined. 

Proposition 4: Let jc = (Jci, ...,x n ) be a tuple of n fuzzy interval-valued variables 
such that xi is any gradual number in the interval defined by gradual numbers \xj ,x~l\. 
Let [zT , z + ] be the smallest fuzzy interval that contains z = f (x \ x„ ) for every x€X, 
where / is the extension of a continuously differentiable real-valued function / :R n i— > 
R. Then 



-/(&) 



and 



max SG3i; ~_/(d)). 



10 



Proof: Finding the minimum function value for f(x\,...,x„) where x £ X = X ,• [x { , xf~ } 
can be stated as a mathematical programming problem: 

Minimize z : f(x) (50) 
subject to c(i) : Jc/— x^ >Q, for2=l,...,« (51) 
c(n + i) : x^—Xj>0, for i = l,...,n, (52) 

where i is a 1 x « vector of fuzzy interval-valued variables. 

Fix a G [0, 1]. Since constraints c(i) and c(n + i) are lower and upper bounds for 
X[ respectively, at most one at a time will be active. We observe that our mathematical 
programming problem satisfies the linear independence constraint qualifications (the 
gradients of the active constraints are linearly independent, which means that the KKT 
conditions are necessary conditions for optimality. Let g(x) = SJ x f(x), and let g(x) be 
the extension of g(x). 

Then for a fixed a G [0, 1], any optimal configuration, 1(a) will satisfy the follow- 
ing: 



(x{a)-Xi[xi(a) -x[(a))-... 




(53) 


...-X„(x„(a) -x~(a)) 




(54) 


+A,„ +1 (I 1 (a)-x+(a)) + ... 




(55) 


... + A, 2n (i;(a)-x+(a)) = 




(56) 


A,,>0, for z = l, 


. . . , 2n 


(57) 


Xi(a) - x[ (a) > 0, i = 1, 


...,« 


(58) 


-jci(a) +I ; + (a) > 0, i = 1, 


...,« 


(59) 


Xi(xi(a)-xf(a))=Q, i=l, 


...,n 


(60) 


X„+i(-Xi(a) +i, + (a)) =0, i=l, 


...,n. 


(61) 



Note that since x (a) = x ! [x i "(a),x / f (a)] is a closed subset of R n , it is a compact 
set. Weierstrass' theorem guarantees that a continuous function attains a maximum and 
minimum on a compact set, so we know f(x) has a minimum on X (a), which has to 

bemk i6 £-/(*(a)). 

The set of points satisfying (l56b-(l6Tb for a fixed a G [0, 1] are exactly the points in 
XjT(a). The minimum of this set will be min- £ ^-f(x(tt)). Taking this result for every 

a G [0, 1] allows us to define z~ by its assignment function SL%- (a) = min ;fG5 ^-/(x(a)), 

which is the same as zT — min~ e ^.- /(f ). 

The proof for z+, is anologous, replacing g(x) in d50l l with — g(x). 

5.2 Example of Fuzzy Interval Analysis via Optimization 

For the sake of illustration, recall our earlier example: 

/(jc)=jc(1-jc). (62) 
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Now let f(x) = x(l —x) be the gradual extension of f(x), and let x be restricted by 
the fuzzy interval X = [xr ,x + ), defined by the following gradual endpoints: 

jT(a) = ia, ae [0,1] (63) 
l+(a) = 1-ia, ae[0,l]. (64) 

Finding the right endpoint of f(x) is equivalent to solving the following optimiza- 
tion problem: 

Minimize x — x 2 (65) 
subject to x — —a > (66) 

-a-l-i>0 (67) 
2 

Now /(x) is continuously differentiable, so we can use Proposition 4 to evaluate 
f(x)onX. 

The Karush-Kuhn-Tucker conditions for ( TTTb are 

(68) 
(69) 

(70) 

(71) 

(72) 

(73) 
(74) 

Now %_j consists of the three triples, (I,A,i,A,2), which satisfy (I68ll-(l74"i>: 

1. {XaMaMa) = (ja,l -a,0) 

2. (xbMbMb) = (1 -ja,0,l -a) 

3. (i c .,X lt .,X 2c .) = (i,0,0). 



1 — 2x — A-i +A,? 











> 


o, j = 


1,2 


1 

i a 

2 

1 - -a -i 
2 

*,i(i--a) 


> 







> 














X 2 {l-^a-x) 









h 


> 


0,A,= 


1,2 



We wish to evaluate / at all three triples to find the minimizer of /. Before pro- 
dding, we note that the first two t 
evaluate feasibility for all a £ [0, 1] 



ceeding, we note that the first two triples in 3y are defined in terms of a, so we must 



1. For (x a ,Xi a ,X2a), all a G [0, 1] are feasible, and f(xb) = jOL— \a 2 ,Va. 

2. For (xb,Xib,X.2b), all a G [0, 1] are feasible, and f{xj,) — ^a 2 ,Va. 



12 



3. In (x e , A-if, A-2 C ), nothing depends on a and/(j) = |,Voc. 

Now (Jc a ,A,ia,X,2a) = (xb,Xu,, A.2&) = (1 — 7«,0, 1 — a) is the minimizer of /, so 
the lower endpoint of the fuzzy interval representing 1(1 — x), is the gradual number 
z~ = 30c- ±a 2 . 

We repeat the steps, minimizing —f(x), to find z + , and find that 3Cyf consists of 
three triples: 

1. {XaMaMa) = (jOC,a- 1,0) 

2. (4,Xi 6 ,X 26 ) = (l-20c,0,a-l) 

3. {xcMcMc) = (3,0,0). 

We wish to evaluate / at all three triples to find the minimizer of /. Before proceeding, 
we note that the first two triples in %_J are defined in terms of a, so we must evaluate 
feasibility for all a <E (0,1]. 

1. For (x a ,X\ a ,X2a), only a = 1 is feasible, and f(x a ) = \ fora= 1 and is undefined 
otherwise. 

2. For (xb, Xib, X2b), only a = 1 is feasible, and f(xb) = \ for a = 1 and is undefined 
otherwise. 

3. In (x c , A-ic, A.2c), nothing depends on a and/(|) = ^, Va. 

Now clearly (x c , X\ c , A-2 C ) = (3, 0,0) is the maximizer of /, so the upper endpoint of 
the fuzzy interval representing x(l — x), is the gradual number tT = \ , Va G [0, 1]. 

Our result, the fuzzy interval [z~,£ + ] = [\& — |oc 2 , ^],Va £ (0, 1] is illustrated in 
Figure |4] 

6 Conclusion 

Viewing a fuzzy interval as an interval of gradual numbers puts at our disposal a num- 
ber of significant tools previously used only for real interval analysis or for analysis 
of a-cuts of fuzzy intervals. We have shown that optimization problems with bound 
constraints can be useful in evaluating non-monotonic functions on fuzzy intervals. 

References 

[1] D. Dubois, H. Fargier, and J. Fortin. A generalized vertex method for computing 
with fuzzy intervals. Proceedings, 2004 IEEE International Conference on Fuzzy 
Systems, pages 541-546, July 2004. 

[2] J. Fortin, D. Dubois, and H. Fargier. Gradual numbers and their application to 
fuzzy interval analysis. IEEE Transactions on Fuzzy Systems, to appear. 



13 



1 



z+ 




alpha 



Figure 4: [z",z + ] = \a\ ±],Voc e [0, 1] 

[3] V. Kreinovich. Probabilities, intervals, what next? optimization problems related 
to extension of interval computations to situations with partial information about 
probabilities. Journal ofGlobabl Optimization, 29:265-280, 2004. 

[4] W. Lodwick. Interval and fuzzy analysis: A unified approach. In Advances in 
Imaging and Electronic Physics, volume 148, pages 75-192. Elsevier, 2007. 

[5] W. A. Lodwick. Constrained interval analysis. Technical Report 138, CCM Report, 
1999. 

[6] R. E. Moore. Interval Arithmetic and automatic error analysis in digital comput- 
ing. PhD thesis, Stanford Univerity, Stanford, CA, 1962. 

[7] J. Nocedal and S. J. Wright. Numerical Optimization. Springer Series in Opera- 
tions Research. Springer, 2 edition, 2006. 

[8] T. Sunaga. Theory of an interval algebra and its application to numerical analysis. 
RAAG Memoirs, 2:547-564, 1958. 

[9] M. Warmus. Calculus of approximations. Bulletin de I'Academie Polonaise de 
Sciences, 3(4):253-259, 1956. 



14 



